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ABSTRACT 

The isotropic distribution of gamma-ray bursts as observed with the 
Burst and Transient Experiment (BATSE) strongly suggests that the bursts 
are at cosmological distances. At such distances, the expansion of the uni- 
verse should redshift the spectra and stretch the temporal structure. Indeed, 
such time dilation has been observed through a variety of analyses and over 
all observed time scales in gamma-ray bursts. We relate the observed peak in- 
tensities, spectral shapes, and time dilation to absolute distance. We include 
the uncertainties in our knowledge of the intrinsic spectrum and correct for 
the coupling between the spectral shape and the temporal structure. Assum- 
ing a go = 1/2 cosmology, the reported time dilation between the dimmest 
BATSE bursts and the bright BATSE bursts (a factor of ~ 2) requires a 
standard candle luminosity of ~ 10 52 erg s _1 , which translates into a red- 
shift of > 6 for the BATSE dimmest bursts rather than a redshift of 1 or 2 
as previously reported. An alternative method to determine the distance to 
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cosmological GRBs is to use the log iV-log P distribution. The large luminos- 
ity and distance determined from the time dilation is inconsistent with the 
observed log TV-log P distribution (which requires a luminosity of 4.6 x 10 50 
erg s _1 and z ~ 0.8) unless there is very strong evolution. Cosmologies with 
qo ^ 1/2 give similar results. The implied distance for the dimmest bursts is 
beyond where galaxies are thought to form. If true, the gamma-ray bursts 
would be orphans: no known objects would have the same distance scale. 
We conclude that either a large fraction (65%) of the observed time dilation 
between the BATSE bright and dimmest bursts is intrinsic to the bursts or 
that there are strong evolutionary effects in the log TV-log P distribution and 
that it is only a coincidence that log TV-log P shows a —3/2 power law at 
high intensities. 

Subject headings: gamma rays: bursts - cosmology: theory 
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1. INTRODUCTION 



Until 1991 there was a strong consensus that gamma-ray bursts (GRBs) 
originated from within our Galaxy. More recently, the Burst and Transient 
Experiment (BATSE) on the Compton Gamma-Ray Observatory discovered 
that there is a dearth of weak events despite the fact that the bursts are 
isotropic (Meegan et al. 1992). This is strong evidence that the bursts orig- 
inate far from the influence of the mass distribution of our Galaxy and even 
the distribution of nearby clusters of galaxies. This has resulted in a great 
debate as to whether gamma-ray bursts are within an extended halo of our 
Galaxy or at cosmological distances. Most observations are either inconclu- 
sive or controversial. The isotropic distribution could be formed by a halo. 
Recently, halos with properties consistent with the BATSE observations have 
been modelled using improved neutron star dynamics (Podsiadlowski, Rees, 
and Ruderman 1994). Harmonically spaced cyclotron lines (Murakami et 
al. 1988, Fenimore et al. 1988) might be a signature of galactic neutron 
stars but it is unclear how such lines can be formed in the super relativistic 
flows required for halo neutron stars (although, see Miller et al. 1991). Even 
though no lines are detected with BATSE, its sensitivity to such lines is not 
sufficient to rule out their existence (Band et al. 1994). A log TV-log P dis- 
tribution gives the number of events (7V) with a peak intensity P (photons 
s _1 cm -2 ) above some value and is sensitive to the distribution of events 
in space. The log N-\og P distribution expected from a cosmological dis- 
tribution has been successfully fit to the combined BATSE - Pioneer Venus 
Orbiter (PVO) observations (Fenimore et al. 1993). However, this is not 
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conclusive since galactic distributions could also fit the data. 

Dimmer cosmological bursts should not only be redshifted compared to 
brighter (presumably nearer) bursts, but the expansion of the universe also 
stretches the bursts in time. Recently, time dilation has been discovered 
in BATSE GRBs (Norris 1994, Norris et al. 1994, Davis et al. 1994). The 
BATSE bursts were assigned to a brightness class depending on their peak 
count rate: the "bright," "dim," and "dimmest" class correspond to 18000 
to 25000, 2400 to 4500, and 1400 to 2400 cts s" 1 in BATSE, respectively. 
(Events with an intermediate count rate are not used since the time dilation 
effects are largest for well-separated classes.) There have been a variety 
of tests that claimed a measured time dilation between bright and dimmest 
events. The total-count test, wavelet-power test, and aligned-peak test find 
a factor of 2 dilation between the dimmest and bright BATSE events (Norris 
1994, Norris et al. 1994). However, Mitrofanov et al. (1994) finds no time 
dilation in BATSE using the aligned-peak test (Mitrofanov et al. 1993) and 
Band (1994) has criticized these tests. Davis et al. (1994) found that a test 
based on the average pulse width required a dilation factor of 1.8, a test 
based on the number of counts within individual peaks required a factor of 
2.2, and a test based on the total number of counts in the bursts required 
a factor of 2.2. Davis et al. 1994 concluded that a time dilation of 2.0 
best explained the data. Fenimore et al. (1995b), using the average auto- 
correlation function, confirmed the time dilation of a factor of 2 between 
the bright and dim BATSE events discovered by Norris et al. (1994) and 
investigated and eliminated possible sources of systematic errors. Recently, 
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Norris et al. (1995) investigated GRB durations and found that the dimmest 
events were dilated by a factor of 2.2 ± 0.2. The dimmest BATSE bursts 
are also softer than the bright bursts as expected for redshifted cosmological 
bursts (Nemiroff et al. 1994). 

The beauty of the time-dilation tests is that they can provide a measure 
of the distance scale to cosmological GRBs that is independent of the log 
TV-log P distribution yet must be consistent with it. The time-dilation tests 
typically use averages from the count rate time history of the bursts as a 
function of brightness, whereas the log TV-log P test uses the number of 
detected events. If these two rather independent tests are consistent, it would 
be strong evidence in favor of the cosmological explanation for the origin of 
GRBs. Indeed, Norris 1994, Norris et al. 1994, and Davis et al. 1994 claimed 
that the observed time dilation of ~ 2 implied that the dimmest bursts were 
at redshifts, z, of ~ 1, consistent with the redshift of the dimmest events 
found from the log TV-log P studies (Fenimore et al. 1993, Wickramasinghe 
et al. 1993). 

In this paper, we find the relationship between time dilation and distance 
correcting for three effects. First, one must assume a spectral shape to 
convert intensity into a distance. In the past, simple spectral shapes (e.g., 
thermal bremsstrahlung or power law) were used and it was assumed that 
one knew the spectral shape at z = 0. In this paper, we will use an average 
over the observed shapes (cf. Band et al. 1993) coupled with their observed 
intensities to determine the distance (redshift) to the events that define the 
spectral shape. Second, the bright class is not necessarily at z = (nor 
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at the redshift of the events that defined the spectral shape) so that the 
time dilation depends on (1 + dimmest) /(I + ^bright) rather than 1 + dimmest 
(Mitrofanov et al. 1994). And, third, since time histories are narrower at 
high energy for a given z, there is a correction based on the redshift of the 
spectrum. 

2. RELATING INTRINSIC LUMINOSITY TO DISTANCE 

As a first step in determining the redshift of a gamma-ray burst from 
observed time dilation, we develop a relationship between standard-candle 
luminosity, Lq, and redshift, z, for different brightness classes. The luminos- 
ity of a source in the detector bandpass depends strongly on the spectrum 
and, since observed spectral shape depends on the distance to the object, the 
intrinsic spectrum of a GRB object must be found. Rather than assume a 
specific spectral shape at z = such as thermal bremsstrahlung, we take as 
our baseline spectra averages over the GRB spectra fit by Band et al. (1993). 
Each such burst has associated with it an observed intensity, Pi^B, and an 
observed spectral shape, (f>i(E), determined by three best-fit parameters pro- 
vided by Band et al. (1993): power law indexes for low and high energies 
and a joining energy. The intensities, Pi t s, m units of photons cm -2 s _1 , 
are given in the BATSE 2B catalog (Meegan et al. 1994) for the range 50 to 
300 keV and are available for 50 of the Band et al. (1993) bursts. 

The observed spectral shape, <f>i(E), will not necessarily come from a 
burst at z ~ 0, especially if Lq is large. Therefore, for a given Lq, Pi b, an d 
4>i(E), we first solve for the redshifts, -z^b, of the events that are associated 
with the spectral shape. Let 4>i(E) be an observed spectral shape for a fixed 
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standard candle luminosity, Lq. The normalization of the spectrum, </>o i, is 
found from: 

£<,= / ^o,^ T-f— hrp— • (!) 

J30 V 1 + z i,Bj 1 + 

The energy range used in calculating the intrinsic luminosity of the source is 
taken as 30 to 2000 keV since we later compare the luminosity to that found 
by Fenimore et al. (1993) from log iV-log P studies. The z = spectral 
shape is <j> ~ ^ , whereas the observed spectrum is 4>(E). The inclusion 

of the 1 + B factor corresponds to the blueshift of the ith baseline burst, 
that is, a burst that we use to define the intrinsic spectral shape. 
The observed intensity of the ith baseline burst is given by: 



, _ Lo C 4>i(E)dE 

i,B " " 



J30 ^ % \l + z i:B J l + zi 

where the comoving distance is 



dE ^Rf z 

B 



(2) 



Ri 



(1 + z %iB )qlH Q 



qoz^ B + (l-q )(l-(l + 2z^ B q )^ 2 ) , (3) 



c is the speed of light, and Hq is the Hubble constant which we take to be 
75 km s _1 Mpc -1 . Except for in section 6, we assume that the deceleration 
parameter, qq, is 1/2, thus, Ri z is 



R ^ = n-L \n (t 1 + z ^ ~ 0- + ^b) 1/2 ) • (4) 

1 + Zi R iSlQ V / 



For a given luminosity, Lq, we numerically determine 1 + for each of 
the Band et al. (1993) baseline bursts using equation (2). The result is 50 
spectral shapes corrected for redshift commensurate with an assumed Lq. 



Norris et al. (1994) used ~ 130 BATSE bursts assigned to the "bright," 
"dim," and "dimmest" classes. Since the publication of Norris et al. 1994, 
another 35 bursts have been accumulated in these classes. For the jth Norris 
et al. (1994) burst in each class, we estimate its redshift, Zj^c, from its peak 
intensity, Pj c f° r each of the Band et al. (1993) spectra (characterized by 

<t>i, z i,B, Lq): 

where P^c is from the BATSE 2B catalog (Meegan et al. 1994) and the 
subscript, C, refers to the brightness class (i.e., dimmest, dim, bright). We 
include (1 + z^c)/ + m this equation to blueshift the observed spectral 
shape back to z = (i.e., divide by 1 + Zi b) and then redshift the spectrum 
to a distance of z^c as seen by us. 

Note that we have parameterized the observed peak intensity in terms of 
the comoving distance. The luminosity distance (Dl = (1 + z)R z ) is appro- 
priate for systems measuring the total bolometric flux whereas GRBs obser- 
vations involve measuring the peak intensity in a finite bandpass. The usual 
bolometric result can be recovered from equation (5). In other applications, 
it is assumed that one knows the intrinsic spectrum, that is, Zi t g = 0. Con- 
verting the number spectrum in the numerator of equation (5) (i.e., 4>(E)), 
into an energy spectrum by multiplying by E and integrating over all energies 
gives the total observed flux: 

L J~E<j>[{l + z)E] dE = L = L 

C E<f>(E) dE AirRl Att(1 + z) 2 R 2 z 4tt£| 



From the 2B catalog we obtained intensity values (Pj c denned over 256 
ms) for 38 of the 45 bursts in the bright class, 37 of the 63 dim bursts, and 
40 of the 59 dimmest bursts. The derived redshifts from equation (5) for 
each of the 38 bright bursts for each of the 50 Band et al. (1993) spectra 
results in 1900 individual redshifts, ^bright? which are averaged. We repeat 
this process for different luminosities and for the BATSE dim and dimmest 
bursts. PVO events make up a fourth class of events even brighter than 
those seen in BATSE. For PVO we use the Band et al. (1993) spectra and 
the typical intensity for a bright PVO event, 200 photons cm -2 s _1 . Figure 
la shows the relationship between standard candle luminosity and redshift 
for the four brightness classes. Also shown on the right hand axis is the 
comoving distance from equation (4). 

The resulting average of 1900 redshifts, zc, accounts for the intrinsic 
spread of spectral shapes, the spread of intensities within a class, and the 
fact that the bursts that define the spectral shape are not at z = 0. For 
example, at a luminosity Lq = 6.3 x 10 51 erg s _1 , the average of the 1900 
redshifts derived from the bright class is ^bright — 1-1? although there are 
individual examples as low as 0.31 and as high as 2.7. The root mean square 
of the distribution of ^bright is 0.41, which we assume is an estimate of 
the intrinsic spread of ^bright- The variation in ^bright comes mostly from 
the variation in Pj,bright, rather than due to the variation in the shape of the 
Band et al. (1993) spectra. Thus, we take as a conservative estimate that the 
number of independent events is 38 (as opposed to 50 or 1900) and the mean 
average ^bright for Lq = 6.3 x 10 51 erg s _1 is therefore found to be 1.1 ± 0.07 
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(i.e., ±0.41/\/38). Similarly, for the same luminosity, Zdim = 4.4 zb 0.18 and 
^dimmest = 5.2 ± 0.20. Note these error bars only reflect the uncertainty in 
the mapping of an assumed luminosity to a redshift; they do not represent 
an uncertainty of a determined value of ^bright , -^dim , or ^dimmest- These 
uncertainties are typical for the curves in Figure la. 

3. ENERGY CORRECTION 
Time dilation due to the expansion of the universe is counterbalanced 
by the narrowing of the time history profile at higher energy: the full width 
half maximum (FWHM) of a pulse as seen in a high-energy bandpass is less 
than that of a low energy bandpass even without any redshift (Fishman et 
al. 1992, Link et al. 1993). Thus, the relationship between time dilation 
and intrinsic luminosity must include a correction factor to account for the 
narrowing of time structure with redshift. Let the dilation factor y account 
for the fact that both the bright bursts are not at z = and the dimmest 
bursts are redshifted at some ^dimmest : 

1 ~l~ ^dimmest / n \ 

V = "1 _l_ y • I 6 ) 

-L ~T~ ^bright 

Notice that y is similar to the spectral correction factor in equation (5) which 
effectively blueshifted the observed spectrum 4>i(E) to z = and redshifted 
the spectrum to a distance 1 + Zj t c- Figure lb shows three different y factors 
each based on a pair of redshifts from Figure 1 function of luminosity 

L . 

The time dilation (Sb-d) °f the dimmest bursts relative to the bright is 
directly related to y. Let W(Ei lt j c ,Ejj t k) be the average FWHM of the time 
structure in the A rth BATSE bandpass (from E^^k to Ejj,k keV) of bright 
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bursts which occurs at an average redshift of ^bright- If that burst occurred 
further away, say at redshift z = ^dimmest, the recorded width in the k th 
bandpass can be derived from the bright bursts using the energy range yE^^ 
to yEu,k- Let W(yEL, k ,yEu, k ) be the average FWHM of a pulse as would 
be seen in the energy range yE^ ik to yEy k keV. The observed time dilation 
of a dimmest burst relative to a bright burst is: 



We refer to W(yEL,k,yEu,k)/W(Ei,,k, Eu,k) as the ^-correction and it 
accounts for the tendency for higher energy peaks to be narrower in individual 
bursts (Link, Epstein, & Priedhorsky 1993). To develop the width function 
dependence on y (i.e., W(yEi,,yE{j)), we fit the observations of average 
pulse widths used by Norris (1994) and Davis et al. (1994) to a model of how 
the pulse width varies with energy. We assume that the time structure can 
be separated from the spectral shape, that is, <f>(t,E) = A(t, E)(f>(E). Let 
R k (E) be the effective area function for the k th BATSE bandpass, then the 
time history observed in the A rth energy band, H k (t), can be modeled as: 



where the energy bands of R k (E) are 25 keV to 57 keV, 57 keV to 115 
keV, 115 keV to 320 keV, and 320 keV to - 1000 keV, corresponding to k = 
1, 2, 3, 4, respectively. Combined channels, such as 25 to 115 keV, use the sum 
of the respective effective area functions (e.g., Ri + 2(E) = R\(E) + R2(E)). 
To construct a profile, we assume a scaling factor, A(t,E), that is based on 



SB-D,k = y 



W(yE L:k ,yE U:k ) 
W(E L , k ,Eu, k ) 



(7) 




(8) 
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an exponential shape: 



A(t,E) = exp 



( 



r(E) 



-t 



(9) 



The spectral shape is not just <f>{E) due to the energy dependence in A(t, E). 
Variations in the spectral shape are not important because they have a mini- 
mal effect on weighting A(t, E) in equation (8). In fact, the variations due to 
A(t,E) are necessary because without spectral evolution through the pulse 
one cannot obtain a narrowing at higher energies. We have found that the 
A(t,E) dependence mimics both the changing intensity within a peak and 
the hard to soft evolution. The decay time of the pulse at a given energy, 
t(E), is assumed to have the form 



where S\ and S2 are unknown parameters and E is in units of keV. The 
motivation for using a power law is twofold. First, the auto-correlation func- 
tion tends to more accurately measure width than the average pulse shape 
and it shows a dependence similar to a power law (Fenimore et al. 1995a). 
Second, the power law fits the data very well (to within tens of milliseconds) 
and agrees well at energies where it can be estimated (see Figure 2). In 
contrast, a linear function for the time decay (i.e., r(E) = S\E + S2) was 
unworkable. In addition, a power law function with a breakpoint does not 
give a significant improvement in the fit. 

We construct a composite profile by averaging the time histories over a 
range of spectra observed by BATSE: 




(10) 
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where each (j>i{E) corresponds to different Band et al. (1993) spectra [n = 



Since the profiles are asymmetric about the peak, we decompose each 
pulse profile into a rise and fall side. The widths (in seconds) are found 
from the profiles as (rise/fall): 0.22/0.44, 0.17/0.32, 0.13/0.27, and 0.08/0.18 
for BATSE channels 1, 2, 3, and 4, respectively. Using a method of least 
squares, we fit the four profile widths generated from equation (11) for both 
sides to these widths and find the four best fit constants for t(E). With the 
constants S 1)T - lse = 0.22, S 1){an = 0.45, S 2 , r i se = -0.40, and S 2 ,faii = -0.39, 
the residual (in seconds) between observed and calculated widths (rise/fall) 
are: 1 x 10~ 5 /9 x 10~ 3 , -7 x 10~ 4 / - 2 x 10~ 2 , 1 x 10- 3 /5 x 10~ 3 , and -9 x 
10 _4 /6 x 10 -3 for channels 1, 2, 3, and 4, respectively. Although no goodness 
of fit can be presented here because we do not know the uncertainties in the 
pulse widths, the calculated widths always match the observed to better than 
a few tens of milliseconds. These residues of the fit are certainly within the 
uncertainties of the widths themselves. 

When the spectrum is redshifted, the photon energy in the equation for 
Hk(t) is shifted by y in those terms that are generated at the burst: 



Notice that the response matrix is unaffected so the only effect on the time 
history is due the to energy redshift in 4>{yE) and A(t, yE). Here, the FWHM 



of H y , k (t) is W(yE L , k ,yE Utk ). Figure 2 shows W{yE L ,yE u )lW{E L ,E u ) 



based on the best fit width functions for BATSE channels 1, 1+2, and 3. 
These width corrections are based on the average pulse shape so, technically, 



54). The FWHM of H k (t) is W(E Ljk ,E Ujk ). 




(12) 
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should only be used to correct time dilation found from the average pulse 
shape test. We believe they can also be used for any test that involves 
linear combinations of peaks (such as the peak-align test) and senses time 
structure the order of a few seconds in a burst. Also shown in Figure 2 
is W(yEj J ,yEjj)/W(Ej J , Ejj) for the auto-correlation test (Fenimore et al. 
1995b) which depends quadratically on the counts. 

One verification that the TF-correction is roughly correct can be seen 
by comparing the ratios of widths constructed from our fitted functions 
at a particular y to the ratio of observed pulse widths for energy bands 
where (yEL,k,yEjj,k) is roughly (Elj, Ejjj). There are two such cases for 
the average pulse energy correction and one for the auto-correlation func- 
tion. For example, from the average fitted pulses, W(Eli, Ejj,\) = 0.65 
sec and W(E Lj2 ,E Uj2 ) = 0.49 sec. Our fitting gives W(2E LA ,2E UA ) = 
0.767W(El,i, Eu,i), which closely matches expectation (0.49/0.65) since a y 
factor of 2 approximately maps the energy range of the first BATSE chan- 
nel (25-57 keV) to the energy range of the second channel (57-115 keV). 
In Figure 2, the solid square represents the observed ratio of W(-Ex i , Ejj\ ) 
to the observed width of W{El^,Ejj^) which does agree with the calcu- 
lated W '(yE L,i,yEu,i) /W (E Eu,i) when y = 2. The solid triangle is the 
observed ratio of W(El£, Ejj, 4) to W(El^, Ejj^), which agrees well with 
W(yE L 2,yEuj)/W(E L j,Eu,3) when y = 2.8. This gives us confidence 
that our width functions are valid at least to y ~ 3. Unfortunately, there 
is not such a point where one can check the width function for the 25 to 
100 keV bandpass that was used in most of the time dilation tests of Norris 

14 



et al. 1994 and Davis et al. 1994. Based on the points where we can check 
them, our width functions seem conservative: the functions actually give a 
smaller correction than if some interpolation of the points were used. It is 
more difficult to establish a width function for the auto-correlation func- 
tion since the auto-correlation function depends quadratically on the counts. 
The open symbol in Figure 3 is a point at which we know the width of the 
auto-correlation function. It, too, is valid to y ~ 3. 

4. RELATING TIME DILATION TO DISTANCE 

Figure la uses the observed intensities to give GRB distances if Lq is 
known. Since we have averaged over spectral shape and burst intensity, 
Figure la only depends on our choice of and Hq. Figure lb is obtained 
directly from Figure la. Figure lc uses the ^-correction from Figure 2 with 
Figure lb to obtain a relationship between the observed time dilation and 
Lq. Thus, the observed time dilation gives Lq which can be used in Figure 
la to give the distance to cosmological GRBs. 

Figure lc is for the BATSE bandpass 25-115 keV. Two different func- 
tions have been plotted corresponding to the bright-dimmest and bright-dim 
y factors from Figure lb. Figure lc assumes the energy-width relationship 
that was used for most of the tests (i.e., 25 to 115 keV) and, technically, 
it can only be used for the average pulse width test. However, given the 
similarity between the various functions in Figure 2, we believe Figure lc 
can be used to estimate the standard candle luminosity for other tests that 
target the same time scale within the GRBs, that is, a few seconds. Figures 
3a and 3b show similar time dilation for the BATSE energy bandpass 25-57 
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keV and 115-320 keV. Figure 3c uses the energy correction from the average 
auto-correlation function (cf. Fig. 2) so provides an accurate mapping of 
time dilation to luminosity for that test. 

Table 1 summarizes the various tests (first column) that have been used 
to detect time dilation. The second column gives the measured time dilation. 
The next two columns ignore the TF-correction: the observed time dilation 
from Table 1 is set equal to y in Figure lb and mapped to an acceptable 
standard candle luminosity using the "Dimmest-Brt" curve. This range is 
reported in column 3 of the table and used with the "Dimmest" curve of 
Figure la to determine the corresponding range of z (column 4). The ±1<t 
ranges for L$o (= Lq x 10 -50 ) are shown in Figure lb for those tests where a 
confidence region has been quoted. The solid horizontal lines at the bottom 
of Figure lb give the acceptable range of the standard candle luminosity 
for the fluence/peak test, the average pulse width test, the average counts 
in a peak test, the total counts test, the duration test, and the average 
auto-correlation test (from bottom to top, respectively). Even without the 
1/F-correction, the observed time dilation indicates a larger standard candle 
luminosity than consistent with the log TV-log P distribution (the dotted 
horizontal line). The best tests (duration and auto-correlation) have an 
average standard luminosity of ~ 20 x 10 50 erg s _1 if the W-correction is 
ignored. In our comparison to the log TV-log P distribution (section 5), 
we will use 2 X 10 51 erg s _1 as an example of the luminosity implied by 
time dilation tests without an ^-correction. Note that the ±lcr range on 
^dimmest (column 4) is usually quite large compared to ^dimmest because of 
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the steepness of Figure lc. This supports our position that one should quote 
time dilation results in terms of the standard candle luminosity rather than 
the z at some ill-defined threshold. 

Energy corrections are clearly required. Columns 5 to 8 in Table 1 incor- 
porate the VF-correction. Some tests (wavelet power, total counts, duration) 
sense time dilation on long time scales within the bursts, and it is less likely 
that Figure 2 is an appropriate ^-correction. These tests cannot be mapped 
to a standard candle luminosity until their energy-width relationship have 
been determined. To include the JF-correction, the observed time dilation 
from Table 1 is set equal to the ordinate of Figure lc (3c for the autocorrela- 
tion test) and mapped to an acceptable range of standard candle luminosity 
using the "Dimmest-Brt" curve. This range is reported in column 6 of Table 
1 and used with the "Dimmest-Brt" curve of Figure lb to obtain the corre- 
sponding y (reported in column 5). The luminosity range is used with the 
"Dimmest" curve of Figure la to obtain the corresponding z and comoving 
distance (i.e., eq. [4]) which are reported in columns 7 and 8 of Table 1. The 
±lcr range for the tests are shown in Figure lc as horizontal solid lines. From 
bottom to top, the lines represent the fluence/peak test, the pulse width test, 
the counts in a peak test, and the auto-correlation test. These tests imply 
that the standard candle luminosity consistent with the time dilation when 
one includes the 1/F-correction is the order of 10 52 erg s _1 . We will use 10 52 
erg s _1 in our comparisons to the log TV-log P distribution (see section 5). 

If the time dilation is due to the expansion of the universe, then it 
must occur equally in all time scales. Indeed, all the tests for which an W- 
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correction can be estimated are consistent. However, the resulting standard 
candle luminosity is much larger than found from the log TV-log P distribu- 
tion. 

5. COMPARISON TO LOG N-LOG P 

One goal of this paper is to investigate the consistency between the dis- 
tance commensurate with the time dilation and the distance commensurate 
with the log TV-log P distribution. There have been a fair number of log 
TV-log P studies, and they have not arrived at consistent results. Most pa- 
pers have concluded that z max is about unity. However, even if one just 
considers results from the same assumed spectral shape, the standard candle 
luminosity concluded by these papers vary by two orders of magnitude from 
< 4 x 10 49 erg s" 1 (Piran 1992) to 4 x 10 51 erg s" 1 (Dermer 1992). An 
additional problem in comparing results arises from using z max since z max 
only has a meaningful value if the threshold it corresponds to is quoted and 
the reader is told how that threshold maps into luminosity. We recommend 
that studies characterize their results in terms of the standard candle lumi- 
nosity in a specified energy range rather than z max . We have used the energy 
range 30 to 2000 keV. Most papers have assumed power law number spectra 
(E~ a ) and the resulting z max and Lq varies strongly with a. For example, 
just varying a from 1.5 to 2.5, Wickramasinghe et al. 1993 found z max could 
vary from 0.5 to 4.0. Associated with these variations with z max were large 
variations in Lq. Given the strong dependence of Lq on the assumed spec- 
tra shape, we suggest that the spectral shape must be accurately modelled. 
Only Dermer 1992 and Fenimore et al. 1993 have used spectral shapes with 
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curvature similar to that observed. However, these studies concluded lumi- 
nosities that were different: Dermer 1992 found Lq = 4 ± 2 x 10 51 erg s _1 
with z max = 1.2 and Fenimore et al. 1993 found L = 5.2+^3 x 10 50 erg 
s _1 in the bandpass 30 to 2000 keV corresponding to z = 0.79 ± 0.05 for 
bursts with peak intensities of 0.7 photons cm -2 s _1 . Dermer 1992 com- 
pared models that ignored threshold effects to the V/V max of the first 126 
bursts observed by BATSE. Threshold effects can only be ignored if the log 
TV-log P distribution is a power law (Band 1993), perhaps contributing to 
the discrepancy. The Dermer 1992 confidence region is large, and it appears 
that it could include the Fenimore et al. 1993 value at a 2a level. We note 
that the Dermer confidence region was based on only 126 events, whereas 
the Fenimore et al. region used bursts from a much larger dynamic range, 
equivalent to 4200 BATSE bursts. 

We will use the Fenimore et al. analysis. This is a crucial decision since 
we will later conclude that the amount of the time dilation is too large to be 
consistent with the distance determined from the log TV-log P distribution. 
It must be recognized that we are not using the log TV-log P study that gave 
the largest distance (Dermer 1992) because of the reasons discussed above. 
We are confident that studies that used only power law spectra (Loredo and 
Wasserman 1994, Mao and Paczyhski 1992, Petrosian, Azzam, and Meegan 
1994, Piran 1992, and Wickramasinghe et al. 1993 as well as our earlier work, 
Fenimore et al. 1992) are undependable. The Fenimore et al. 1993 log TV-log 
P study combined events from BATSE and PVO. Due to its small size, PVO 
only saw the brightest events, but PVO had a longer lifetime (~ 14 yr), larger 
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on-time (90%), and a larger field of view (4-7T since PVO had very little Venus 
blockage). It will take BATSE ~ 30 more years to observe as many bright 
bursts as PVO. We apply a uniform selection criterion for the included events 
and accommodate the differences in instrumental response in calculating 
intensities. The threshold effects are accounted for. In Fenimore et al. 1993, 
we assumed that the intrinsic spectral shape is (f>(E) = E~ l exp(— P/350) 
such that 

/.2000 

L = E<j)Q<j)(E)dE (13) 

^30 

and related z to intensity by 

p L C m + z j£ )E]dE 

/ 30 E<f>(E)dE ^R jz 

(Note in Fenimore et al. 1993 we used an energy range that was common to 
both BATSE and PVO, 100 to 500 keV.) We fit the observed number of 
bursts in a range of peak intensities to 

A7V(P a to P 2 ) = / — — dV . (15) 

Here, p is the rate-density of events per comoving volume (bursts yr _1 
Gpc -3 ). In contrast to steady sources (such as galaxies), equation (15) re- 
quires the extra 1 + z in the denominator because the time between bursts is 
dilated with the expansion of the universe. For go = 1/2, dV has the simple 
form 

dV = 4vrr 2 dr. (16) 

Such an analysis effectively assumes that we know the intrinsic spectral 
shape (i.e., at z = 0) when, in fact, we have only observed it for the bright 
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bursts which would be at a substantial z if Lq is large. Thus, in this pa- 
per we make two improvements to our previous log TV-log P study. (These 
improvements have not been utilized in any other study that we are aware 
of.) First, when we do the fitting using equation (15) we average over a set 
of spectral shapes that span the parameter space for GRBs, that is, we use 
the 50 Band et al. 1993 spectra in equations (1) and (5). Second, we intro- 
duce a new parameter, zb, the redshift of the events that define the spectra 
shape: compare equation (1) with (13) and (5) with (14). We use an addi- 
tional relationship (eq. 2) with additional information (P^b) to determine 
Z B- (When using eq. [5] for the log TV-log P studies, we use an integration 
range of 100 to 500 keV in the numerator; and when we use eq. [5] for Fig. 1 
or 3, we use a range of 50 to 300 keV.) 

The histogram in Figure 4 shows the log TV-log P distribution of the 
combined, uniformly selected events from PVO and BATSE. Table 2 provides 
the intervals (Pi and P2) and the number of observed events (ATV &<.) used 
for the fitting of the differential log TV-log P to equation (15). The fitting 
used a % 2 statistic defined as 



where A s normalizes the observed number of events in a photon s cm 
range to events yr" 1 (47T)" 1 . For PVO, A a = 0.09072 and for BATSE, 
A s = 4.237 (see Fenimore et al. 1993). We did not use data below where 
either instrument has threshold effects. Note that this uses the model rather 
than the observations to estimate the variance on the observations. Table 
2 can be used to fit any model to the combined BATSE-PVO log TV-log 




{A. s AN obs - ANjP, toP 2 )) 2 
^ATV(Pi to P 2 ) 



(17) 
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P. The results of the re-analysis is Lq = 4.6^Q'gg X 10 50 erg s 1 in the 
bandpass 30 to 2000 keV with p = 24 Gpc -3 yr _1 , rather than 5.2 x 10 50 
erg s _1 as reported in Fenimore et al. 1993. The solid line is the expected 
distribution given a best-fit standard candle luminosity of 4.6 x 10 50 erg s _1 
the corresponding % 2 is 9.1 with 9 degrees of freedom. The dashed line in 
Figure 4 represents the expected log TV-log P distribution of events given a 
standard candle luminosity of 10 52 erg s _1 , the value consistent with the time 
dilation tests; it poorly matches the observed distribution from BATSE and 
PVO (x 2 is 510). The dotted line is for L = 2 x 10 51 erg s , the standard 
candle luminosity obtained for time dilation without any l/F-correction. It, 
too, poorly matches the observed BATSE-PVO log TV-log P distribution (% 2 
is 60). 

Perhaps evolution can explain the discrepancy (Norris et al. 1995). The 
expected log TV-log P distribution can be modified by either number evo- 
lution (i.e., p = p(z)) or luminosity evolution (i.e., L = L(z)). Spectral 
evolution has little effect on the expected distribution. We use the func- 
tional shapes L = Lq(1 + z) PL and p = po(l + z) p " which were assumed 

3 /2 

by Fenimore et al. 1993. However, since TV scales as pL for most of the 
parameter range, we effectively fit with (1 + z) ppJ <~i pL . In Figure 5a, we show 
the evolution necessary (p p + to achieve agreement with the BATSE- 

PVO log TV-log P as a function of standard candle luminosity. The fits are 
unacceptable for p p + t^TPl greater than ~ 2.4. Of course, other functional 
forms for the evolution might be able to provide acceptable fits. The required 
evolution is quite strong and a successful cosmological GRB model must ac- 
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commodate it. For example, assume p p + \vl = 2.4 even though it provides 
insufficient evolution to obtain a best fit (cf. Fig. 5a). The time dilation re- 
quires z to be at least 6 (cf. Table 1). Thus, if one has density evolution, the 
density of progenitors must be a factor of ~ 100 higher at the distance of the 
dimmest bursts. It requires a factor of ~ 10 higher density at the distance of 
the bright BATSE bursts. In the context of a colliding neutron star model 
for GRBs, this would require that neutron star production was much higher 
in the distant past (which is not inconceivable). Alternatively, if luminosity 
evolution occurs, pi is 1.6 and the objects must be ~ 25 times brighter at the 
distance of the dimmest events. Since neutron stars have a relatively small 
range of masses, it seems unlikely that colliding neutron stars would release 
more energy just because the universe is younger. However, the process that 
converts the released energy to gamma rays is uncertain and could depend on 
epoch. Epstein et al. (1993) suggested that GRBs are upscattered ambient 
photons from colliding stars near an AGN. AGNs are brighter in the distant 
past so, perhaps, the reason GRBs have luminosity evolution is that there 
are more ambient photons in the vicinity of the colliding neutron stars. 

Perhaps more troublesome is that the evolution required for consistency 
between the time dilation and the log TV-log P distribution coincident ally 
produces a —3/2 power law in the PVO log TV-log P distribution. Figure 5b 
gives the % 2 f° r fitting just the 6 PVO bins above 20 photons cm 2 s _1 and the 
BATSE bin above 7.08 photon cm 2 s _1 (see Table 2) with equation (15). This 
is the region of the log TV-log P distribution that follows a -3/2 power law. 
Low luminosities give an acceptable % 2 of 5.0 (with 5 degrees of freedom). 
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The x 2 grows quickly for luminosities greater than 10 51 erg s _1 . For the 
luminosity consistent with time dilation (with the TV-correction), the % 2 has 
grown to 30. Evolution is required to obtain an acceptable fit in the -3/2 
portion of the distribution. Thus, the shape of the dashed curve in Figure 4 
is inconsistent with the —3/2 portion of the log iV-log P distribution. Strong 
evolution can modify the shape to agree with the -3/2 portion but then it 
must be considered a coincidence that the log TV-log P distribution has a 
portion that appears to be from a homogeneous distribution. The chances 
that such a coincidence would occur is roughly the probability of obtaining 
a x 2 °f 30 with 5 degrees of freedom. 

In previous studies it was usually assumed that we knew the intrinsic 
shape of the spectrum (i.e., eqs. [13] and [15] were used rather than eqs. [1] 
and [5]). The inclusion of zb improves the comparison between the time 
dilation and the log TV-log P distribution. If zb is assumed to be zero, then 
Figure 1 would actually give even a larger Lq and ^dimmest- The inclusion of 
the zb terms moderates variations in Lq. Even though ^dimmest is quite large 
(> 6), the corresponding Lq does not vary much from a low z value. The 
reason is that (1 + 2dimmest)/(l + -^bright) is limited by the dynamic range of 
the observed GRBs peak intensities, which is small (about a factor of 30 for 
BATSE). 

6. OTHER COSMOLOGIES 

The luminosity implied by the time dilation has many fewer faint events 
than observed (cf. Fig. 4). Strong evolution will allow some reconciliation 
of the log TV-log P distribution and the time dilation. Another possible way 

24 



to reconcile the difference is to assume different cosmologies. Cosmologies 
with low qo will produce relatively more faint events for the same standard 
candle luminosity. However, we will show that small go requires a much 
larger luminosity such that the time dilation is not consistent with the log 
iV-log P distribution. For qo near zero, equation (3) becomes 

cz . z, , 

R " = s^' 1 + 5»- (18 » 

Using equation (18) in equations (2) and (5), we produce Figure 6 which 
is similar to Figure 1. The redshift, z, varies more slowly with luminosity 
than in the case of qo = 1/2. At a given luminosity (say, 10 52 erg s _1 ), 
The z of the dimmest events is not much different than the z of the bright 
events. Thus, y is relatively smaller for qo = compared to go = 1/2. If 
one ignores the W-correction and identifies the observed time dilation of 2.0 
with y in Figure 6b, then the required luminosity is ~ 2.8 x 10 52 erg s _1 . 
When qo ^ 1/2, dV in equation (15) takes the form 

\H J q*(\ + 2)3(1 + 2<, 2)l/2 

(see eq. [2.56] in Kolb & Turner 1992 ). Using equation (19) in equation (15) 
gives a % 2 of 314 for Lq = 2.8 X 10 52 erg s _1 . Since y is always small (due 
to the relatively small differences in z in Fig. 6a), there are no reasonable 
luminosities that produce a time dilation as large as observed (Fig. 6c). 

The other extreme case is to assume that qo = 1.0. Although unlikely, 
it does allow better agreement between the log iV-log P distribution and the 
time dilation. In this case, equation (3) becomes 

= J^ttty <20) 
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For qo = 1, we find that the log iV-log P distribution is best fit by a luminosity 
of 3.7 x 10 50 erg s _1 , with x 2 = 9.1 for 9 degrees of freedom. Using equation 
(20) in equations (2) and (5), one can repeat the analysis of Figure 1 and 
determine the standard candle luminosity that is consistent with the observed 
time dilation. One finds for qo = 1 that a luminosity of 7 X 10 50 erg s _1 has 
a time dilation of 2.0 with no PF-correction and a value of 2.75 x 10 51 erg s _1 
including the ^-correction. The assumption of go = 1 allows the dimmest 
BATSE events to be at a lower z, about 5. Figure 7 is similar to Figure 4 
except it assumes qo = 1. A luminosity of 7 x 10 50 erg s _1 fits the log TV-log 
P distribution with a % 2 of 18 and a luminosity of 2.75 X 10 51 erg s _1 fits 
with a x 2 °f 180. Thus, the luminosities that fit the time dilation are closer 
to that which fits the log TV-log P distribution but are still inconsistent. 

In Figure 5, the dotted curves are for go = 1- At 2.75 x 10 51 erg s _1 , 
substantial evolution is required (p p + |pl = 3.2, Fig. 5a). If one just fits to 
the -3/2 portion of the log TV-log P distribution, the luminosity consistent 
with the time dilation fits with a % 2 of ~ 12 (see Fig. 5b). Since the best 
fit x 2 f° r the "3/2 portion is 5 and there are two fit parameters, the time 
dilation luminosity is approximately at the 2cr contour of acceptable log TV- 
log P luminosities. Thus, one can reconcile the time dilation and log TV-log 
P distribution if qo = 1 and there is substantial evolution (p p + §Pl = 3.2). 
In that case, it is a 2a coincidence that we see a -3/2 power law at high 
intensities. 

7. DISCUSSION 

All tests show time dilation of ~ 2 between the BATSE bright bursts 
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and the dimmer bursts (Norris 1994, Norris et al. 1994, Davis et al. 1994, 
Fenimore et al. 1995b, although see Mitrofanov et al. 1994). The key question 
is what distance does a time dilation of 2 imply and is it consistent with the 
log TV-log P distribiution. Norris et al. (1995) estimates Zdimmest based on a 
consistency argument. Norris et al. (1995) uses Sb-d = 2 ~ y as a first order 
approximation to estimate the ^-correction and obtains a ^-correction for 
the average pulse width test equal to 0.85 which corrects y to 2.35. The Norris 
et al. (1995) conversion from y to z uses equation (6) and assumes zs r t = 0.3 
(the log TV-log P value for the bright BATSE events) and concluded that 
Zdtmmest ~ 2. This is close to being self-consistent, for z dtmmest = 2,z hrt 
is actually ~ 0.6. But, as z^ rt get larger, the PF-correction and Zdimmest 
get larger causing Zbrt to be larger. Thus, Norris et al. (1995) obtained a 
low value for Zdimmest primarily because z^rt was assumed to be 0.3. The 
purpose of this paper is to formulate the relationship between time dilation 
and distance including all known effects. Rather than an iterative process, 
we calculates the function in Figure lc which maps Sb-d directly to Lq and, 
therefore, Zdimmest and y. A Sb-d of 2 gives a y of 3 and the ^-correction is 
~ 0.68. We use a different conversion from y to z then did Norris et al. (1995). 
For example, our conversion recognizes that the intensity and spectral shape 
used to define z vs. Lq were not from z = (note the Zf, terms in the eqs. [1] 
to [5]). Thus, even if one accepts the Norris et al. (1995) y value of 2.35, the 
corresponding Zdimmest found by Figure 1 is 4, not 2. When one uses Figures 
lc and la, the Zdimmest corresponding to, for example, Sb-d = 2lg'25 ^ s 
5.5^2 4 ( see Table 1). Thus, for the same observation, different analysis 
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methods have resulting dimmest ' s that vary from 2 to 5.5. At a low value 
of z = 2, the log TV-log P distribution result (z = 0.8, Fenimore et al. 1993) 
might be within the statistics or modest evolution could produce consistency 
(Norris et al. 1995). However, this distance scale for the dimmest events 
was arrived at by assuming that the bright events were at the z consistent 
with the log TV-log P distribution. Although reasonable if z\, rt is small, a 
more detailed analysis is required if ^ mmes t^0.5. At z = 5.5, even extreme 
evolution has difficulty obtaining consistency (see section 5). 

Even if the ^-correction term is ignored, time dilation of ~ 2 implies 
^dimmest is ~ 2.5 or a standard candle luminosity of 2 x 10 51 erg s _1 . This 
result depends only on the conversion between peak intensity and standard 
candle luminosity (i.e., Fig. la). To be consistent with the log TV-log P distri- 
bution, such a high luminosity requires some evolutionary effects. When the 
required VF-correction is included, an observed time dilation of ~ 2 actually 
maps into a ^dimmest of >6 that is, a large standard candle luminosity (10 52 
erg s _1 ) requiring substantial evolution to be consistent with the log TV-log P 
distribution. This result depends on our estimate of the TF-correction term 
in equation (7), which is shown in Figure 2. A time dilation of 2 implies that 
y is about 2.9. The ^-correction term is reasonably well determined to at 
least y equal to 3 (see solid points in Fig. 2). Thus, we do not think that the 
TF-correction term is falsely implying a large Lq. When y ~ 2.9, the average 
redshift (z) of the BATSE bright bursts is ~ 1.2 and the average redshift for 
the dim events is ~ 5.6 (cf. Fig. la). Thus, when there is net time dilation 
of ~ 2, the time dilation due to the expansion of the universe is responsible 
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for a factor of 2.9, but that is counteracted by the narrowing of the pulse 
profiles with energy which contributes a factor of 0.68. 

Beaming can reduce the overall energy requirements to be less than 
10 52 erg s _1 . However, beaming does not change our estimate of ^dimmest; 
it only scales the reported luminosity by fl/Air where fl is the angular size 
of the beam. Beaming does not affect our comparison to the log TV-log P 
distribution, since Lq scales in the same way in both analyses. 

For many models, a z of > 6 is uncomfortably large. For example, to be 
consistent with the neutron star merger model (Meszaros and Rees 1992), one 
must assume a very early phase with a much higher density of neutron star 
pairs. Since galaxies are thought to form at about z ~ 3, it is not clear why 
there should be so many neutron star pairs so early On the other hand, the 
greater distance would help resolve the "no host" problem. Schaefer (1992) 
pointed out that small error boxes for bright (nearby) GRBs are devoid of 
bright extragalactic objects such that M31-like galaxies could not be the host 
galaxy. Fenimore et al. 1993 showed that it was difficult to reconcile the log 
TV-log P distribution with any host galaxy if GRBs followed the luminous 
mass. Certainly bright AGNs were ruled out. If the greater distance implied 
by the observed time dilation is correct, then it is not surprising that we 
do not see the host galaxies. The no-host problem does not exist if GRBs 
are associated with galactic centers since the number of galaxies (in contrast 
to the luminous mass) could be very large for small, virtually undetectable, 
galaxies. 

There are three possible explanations that could reconcile the excessive 
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time dilation with the log iV-log P distribution. First, the time dilation 
could, indeed, be a result of the expansion of the universe and the bursts 
come from extreme distances (z beyond 6). The lack of normal galaxies 
within the error boxes would no longer be a problem, although the bursts 
would be orphans: no know objects with the same distance scale. Either 
strong density evolution or luminosity evolution would be required. The 
luminosity evolution might result from models that utilize the ambient radi- 
ation of AGNs (cf. Epstein et al. 1993). Under this explanation, the observed 
-3/2 power law in the log TV-log P distribution must be considered a coin- 
cidence rather than indicating that some of the bursts are from a nearby, 
homogeneous population. 

Second, the log TV-log P distribution could indicate a valid distance 
scale, but, therefore a large fraction of the observed time dilation is intrinsic 
to the bursts. We tend to favor this explanation since it is probably easy to 
have intrinsic time stretching of the bursts. The distance indicated by the 
log TV-log P distribution implies that the dimmest GRBs should only have a 
time stretch of ~ 1.3 (cf. Fig. lc) rather than the observed 2. Thus, 65% of 
the time stretching must be intrinsic in order to be consistent with the log 
TV-log P distribution. Under this explanation, it would be difficult to use 
the time dilation as an argument that GRBs are cosmological since perhaps 
all of the time dilation could be intrinsic. 

Indeed, Band (1994) and Wijers & Paczyhski (1994) have argued that 
such time dilation could easily be intrinsic to GRBs. Clearly, the most 
suspect assumption that Norris et al. (1994) and we have made is that peak 
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intensity indicates distance. Just based on the variation of peaks within 
individual bursts, it is reasonable that GRB peak intensities have a very 
broad distribution even for sources that are at the same distance. If bright 
bursts intrinsically have shorter time scales, then we could just be seeing the 
result of a luminosity function. This could arise easily if, for example, fluence 
(~ luminosity times duration) was standard rather than peak intensity, see 
Band (1994) and Wijers & Paczyhski (1994). All reports of time dilation have 
assumed a standard candle luminosity and, if bursts have, in fact, a standard 
candle fluence, then the reports of time dilation are probably artifacts of 
the standard candle luminosity assumption. Unfortunately, total fluence is 
harder to measure and the trigger efficiency for fluence is so poorly known 
that a standard candle fluence analysis is probably not possible for time 
dilation. 

Third, all of the time dilation tests that we have used (Norris 1994, 
Norris et al. 1994, Davis et al. 1994, Norris et al. 1995, Fenimore et al. 
1995b) are based on the same set of events. Although all efforts to identify 
systematic effects have been taken and statistical variations accounted for, 
perhaps there is an unknown systematic effect that gives the appearance of 
a correlation of intensity and time scale. Not all bursts were used, some had 
data gaps that made them unusable and these tend to be the longer bursts. 
The shortest events (less than 1.5 s duration) were omitted, although it is not 
clear why that would induce a correlation. It is our judgment that systematic 
effects have been accounted for. 

Finally, we make the following recommendations concerning time dila- 
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tion tests. First, one must use yW(yEL, yEjj) jW (El, Ejj) to properly relate 
the observations to a distance scale rather than 1 + ^dimmest- Although the 
differences are small they can cause a large effect. Second, for each test, 
the bright events must be used to provide W(El, Ejj). Third, ^dimmest ^ s 
rather nebulous because it depends on the threshold of the instrument. One 
should quote the equivalent Lq in a specified energy range. Error bars on 
^dimmest often are large because z diverges at distances near 2c/ Hq. Fourth, 
the assumption of simple spectral shapes, such as power laws or thermal 
bremsstrahlung, is unwarranted. One should use averages over the Band et 
al. 1993 spectra. And, fifth, one cannot assume that the events that define 
the spectral shape that are used in either the log TV-log P studies or the con- 
version of time dilation to distance are at z = 0. One needs to use equations 
(1), (2), and (5) to properly relate peak intensity to distance. 
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FIGURE CAPTIONS 

Fig. 1: a) Average 1 + z and comoving distance versus standard candle 
luminosity Lq for various brightness classes: the dotted, dashed, and dot- 
dash curves are for the bright, dim, and dimmest classes of Norris et al. 
1994, respectively. The long dash curve is for typical bursts from PVO. 

b) The ratio of the 1 + z factors for three pairs of brightness classes: the 
dotted, dot-dash, and long dash curves are dimmest /bight, dim/bright, and 
bright/PVO, respectively. If there is no ^-correction, Fig. lb gives the 
relationship between observed time dilation and standard candle luminosity. 
The solid horizontal lines represent acceptable ranges (±ler) from various 
time dilation tests assuming no 1/F-correction. The horizontal dotted line at 
4.6 x 10 50 erg s _1 represents the acceptable range of the standard candle 
luminosity consistent with the log TV-log P distribution. 

c) Observed time dilation versus luminosity including the ^-correction for 
the bandpass 25-115 keV. The solid horizontal lines represents the acceptable 
ranges of luminosity (±ler) from various time dilation tests if one includes 
the ^-correction. If GRBs have the luminosity associated with the log TV-log 
P distribution, then one would expect to see only a factor of 1.3 time dilation 
between the dimmest and brightest BATSE GRBs rather than a factor of 2 
(cf. panel c). 

Fig. 2: The W-correction function, W(yE L ,yEu)/W(E L , E v ), as a 
function of y as determined by fitting a model of the time variation to the 
average pulse shapes. The solid, dotted, and dashed curves are for the 25 to 
115 keV, 115 to 320 keV, 25 to 57 keV bandpass in BATSE for the average 
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pulse shape test. The long dashed curve is for the average auto-correlation 
test in the BATSE 115 to 320 keV bandpass. The points are values where the 
function can be estimated directly from the observations. The filled square 
is for the 25 to 57 keV bandpass at y = 2, the filled triangle is for the 115 to 
320 keV bandpass at y = 2.8, and the open triangle is for the 115 to 320 keV 
bandpass at y = 3 for the average auto-correlation correction factor. These 
points demonstrate that our model is valid to at least y ~ 3. 

Fig. 3: a) Time dilation versus luminosity for the 25 to 57 keV bandpass 
of BATSE. The curves are labelled as in Fig. lc. 

b) Time dilation versus luminosity for the BATSE bandpass 115-320 keV. 
The solid horizontal line is the acceptable range of luminosity (±1<t) from 
the average auto-correlation test for time dilation assuming no PF-correction. 
The dotted horizontal line at 4.6 x 10 50 erg s _1 represents the acceptable 
range of the standard candle luminosity consistent with the log TV-log P 
distribution. 

c) Time dilation versus luminosity suitable for the average auto-correlation 
test. The Dimmest-Brt and Dim-Brt curves use the > 320 keV energy range 
of BATSE. The PVO-Brt curve uses the > 115 keV energy range of BATSE 
and the 100 - 2000 keV energy range of PVO. The solid horizontal line is 
the acceptable range (±ler) from the average auto-correlation test for time 
dilation including the 1/F-correction. 

Fig. 4: The log TV-log P distributions for various standard candle 
luminosities compared to the BATSE and PVO observations assuming go = 
1/2. The solid curve (L = 4.6 x 10 50 erg s -1 ) provides the best fit to the log 
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iV-log P observations. Only data above the threshold for the instruments 
were used: 1.0 photons s _1 cm -2 for BATSE and 20 photons s _1 cm -2 for 
PVO. The dotted curve (Lq = 2 x 10 51 erg s _1 ) corresponds to the luminosity 
expected for an observed time dilation of 2.0 without an VF-correction. The 
dashed curve (Lq = 10 52 erg s _1 ) corresponds to the luminosity expected 
for an observed time dilation of 2.0 with an ^-correction. Although strong 
evolution might be able to force the distributions with large luminosities into 
agreement with the observations, it implies either a higher density of GRB 
progenitors and/or more luminous events in the distant past and that the 
observed —3/2 power law in PVO is a coincidence. 

Fig. 5: a) The density and/or luminosity evolution necessary to achieve 
agreement with the BATSE-PVO log TV-log P distribution as a function of 
the standard candle luminosity. The density evolution is assumed to have 
the form (1 + z) Pp and the luminosity evolution has the form (1 + z) PL . 
The extent of the curve represents the range over which acceptable fits are 
possible. The solid line is for cosmologies with = 1/2 and the dotted line 
corresponds to qa = 1. 

b) The x 2 for fitting the -3/2 power law portion of the BATSE-PVO log 
TV-log P distribution to cosmological models with different standard candle 
luminosities. For luminosities greater than ~ 2 x 10 51 erg s _1 , it should be 
considered a coincident that the log TV-log P distribution shows a -3/2 power 
law since it arises because of the evolution, not because of a homogeneous 
distribution of sources. 

Fig. 6:The average 1 + z, y, and time stretching as a function of the 
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standard candle luminosity, Lq, if = 0. The labelling is the same as in 
Figure 1. When go = 0, the bright BATSE events have values of z that 
are closer to the those of the dimmest events (see panel a). Thus, the time 
stretching without VF-correction (panel b) is smaller than if go = 1/2- The 
time stretching with 1/F-correction (panel c) is smaller yet such that no lu- 
minosity within the range of our assumptions can give a time stretching as 
large as observed (i.e, 2.0). The horizontal dotted line at 4.6 x 10 50 erg s _1 
represents the acceptable range of the standard candle luminosity consistent 
with the log TV-log P distribution. 

Fig. 7: The log TV-log P distributions for various standard candle 
luminosities compared to the BATSE and PVO observations assuming qo = 
1. The solid curve (Lq = 3.7 x 10 50 erg s _1 ) provides the best fit to the log 
TV-log P observations. The dotted curve (Lq = 7 x 10 51 erg s _1 ) corresponds 
to the luminosity expected for an observed time dilation of 2.0 without an 
1/F-correction. The dashed curve (Lq = 2.75 X 10 51 erg s _1 ) corresponds 
to the luminosity expected for an observed time dilation of 2.0 with an W- 
correction. 
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TABLE 1 

Distances to Cosmological GRBs from Time Dilation Tests 







No W(E L ,Eu) 




with W(E L , 




Test 


S 


L50 


^dimmest 




L50 


^dimmest 


d (Gpc//i 75 ) 


Peak Align 1 


~ 2.25 


22 


2.5 


3.6 


126 


9.5 


5.5 


Fluence/Peak 1 


9 n+0.50 
z - u -0.25 


13i 22 


l-°-o.e 


3.0 






4.9± ;° 8 


Wavelet 1 


- 2.25 


22 


2.5 










Pulse Width 2 


1 o+0.65 




1 4+1-8 


2.5 


36l 118 


c + 9-4 
o.o_ 2 _7 


4.2i;;5 


Cts in Peak 2 


9 9+0.72 
z - z -0.44 




2.4±2;f 


3.5 


H3153 3 


8.6l 12 4 3 


5.4+^ 


Total Cts 2 


9 n +0.62 
z - u -0.47 


1 o + 28 
i,3 -10 












Duration 3 


2.2 ± 0.2 


20t 9 7 


9 4 +o.6 
z - 4 -0.6 










Auto-Corr 4 


2.0 ± 0.2 


2li 12 


9 c + 0.8 
z -°-0.7 


3.0 


120lg 


9-o±!:l 


c c+0.4 
°-°-0.7 



Morris et al. 1994 
2 Davis et al. 1994 
3 Norris et al. 1995 
4 Fenimore et al. 1995b 
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TABLE 2 

Data for Fitting to the BATSE-PVQ log TV-log P 



Pi 


P 2 


Satellite 


AN obs 


1.00 


1.26 


BATSE 


18 


1.26 


1.78 


BATSE 


20 


1.78 


2.82 


BATSE 


12 


2.82 


7.08 


BATSE 


26 


7.08 


31.6 


BATSE 


10 


20.0 


25.1 


PVO 


41 


25.1 


31.6 


PVO 


20 


31.6 


39.8 


PVO 


27 


39.8 


56.2 


PVO 


21 


56.2 


100 


PVO 


23 


100 


1000 


PVO 


14 
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